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OF A TVro-MASS SYSTEM WITH AN APPLICATION 
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STMIAEY 


Solutions of Impact of a rigid prismatic float connected ty. 
a massless spring to a rigid upper mass are presented. The 
solutions are "based on hydrodynamic theory which has been experi- 
mentally confirmed for a rigid structure. 

Equations are given for defining the spring constant and 
the ratio of the sprung mass to the lower mass so that the two- 
mass system provijies representation of .the fundamental mode of, 
an airpl^e wing. The forces calciilated are more accurate than, 
the forces which would be predicted for a rigid alrfratiie since 
the effect of the fundamental mode on the hydrodynamic force is. 
taken into account. The response of the two-mass system gives , 
the response of the. represented mode and, although no provision 
is made for taking into account the effect of secondary modes on 
the hydrodynamic force, means are indicated whereby the resvilts 
may be used to approximate the response of modes other than the 
fundamental mode. 

Time histories of the hydrodynamic force and structural 
response are given for wide ranges of mass distribution and ratio 
of natural period to the period of the impact. By use of non- 
dimensional coefficients these results are made applicable to 
different combinations of velocity, weight, angle of dead rise, 
and fluid density. Although the equations permit solutions for 
different combinations of flight-path angle and trim, an approxima 
tion is given for correcting' the results for the combination for 
which solutions are given to other conditions within a narrow 
range indicated. to. be of • primary Interest to the design engineer. 

In a comparison of the theoretical data with data for a 
severe flight-test landing impact, the effect of the fundamental 
mode on the hydrodynamic force is considered and response data 
are compared with experimental data. Consideration of the 
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fundamental mode alone fails to account for the fact that during the 
impact partial failure of the inboard -engine moxmts occurred but 
use of the theoretical solutions to approximate the effects of 
further ving torsion leads to substantial agreement. 


INTRODUCTION 


In recent years the development of large airplanes has caused 
the elastic behavior of airframe structures during landing impact 
to become important. The work which has been done on this problem 
has been handicapped by lack of proper knowledge of the time 
history of applied ground reaction. This situation has been 
particularly acute for seaplanes because of difficulties in 
measuring the hydrodynamic force, the^ seaway, and the manner of 
contact with the seavray. . 

In order to facilitate the interpretation of flight data and 
to lead to the prediction of design loads on a rational basis, a 
theoretical hydrodynamic .study was made and tests of a rigid float 
were conducted at the Langley impact basin* Since the results of 
these tests agree with the theoretical results for wide .ranges of 
the pertinent variables in numerous force time histories (refer- 
ence 1), it is assumed that the theory may also be used in con- 
sidering the effect of the upper- structure elasticity of a seaplane 
on the motion and force characteristics of the hull proper, which 
is assumed to be rigid. 

The bending of wings during impact, which for modem flying 
boats is the primary structural action, is considered in the 
present paper by reducing the fundamental mode to an equivalent 
two-mass system. The results are presented in a form suited to 
general application. and are compared with experimental results 
for a partlc\ilar case. The equations showing the method of 
solution are included in appendix A and a sample data sheet is 
given as table I . 


SYMBOLS 


tjj time required for one -fourth cycle of natural vibration 

tj^ time between initial contact and maximum hydrodynamic 

force for rigid structure 

t ' time elapsed after Initial contact 
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lower, or hiill, mass of two-mass system 
m0 sprung mass of two-mass system 

m gross mass (w/g or mg + raj,) 

¥ gross wel^t 

g acceleration of gravity 

K spring constant of spring connecting mg and force 
. per unit dsf lection 

nj acceleration normal to water surface of nodal point of 
elastic .system, multiples of acceleration of gravity; 
for two-mass system, acceleration of center of gravity 

' r-o oscillatory acceleration of hull about center of gravity 
-of two-mass system or nodal point of represented mode, 
multiples of the acceleration of gravity 

V ■ resultant velocity at instant of contact wi^ih water surface 

p mass density of fl\ild 

T angle of trim,- angle of hull keel with respect to plane of 
water surface 

r flight-path angle; angle between fli^t path and plane of 
water surface 


0 angle of dead rise 

C.J. nondimensional time coefficient 

C7 nondimensional load-factor coefficient 


nondimensional draft coefficient 




^%ax Instant of maximum acceleration 


fj^ natural bending frequency 

Where units are not given, any consistent system of unite may be 
used . • 
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The hydrodynamic theory used in the present paper is the same 
as that developed in references 1 and 2. A haslc differential 
equation vhich gives the instantaneous force in terms of the 
instantaneous position and motion of the float, is given in refer- 
ence 2. This equation in the present paper is used to determine 
the effect of airframe elasticity in altering the motion ^d force 
time history (appendix A) . The solution is heieed on the assumption 
that the float does not change trim during impact. In this connec- 
tion the pitching moment may he large, hut the time of the 
Impact is short enough. to warrant (at the present steige) neglect 
of the resu3.tlng angular velocities and displacements. 

The solution presented herein is for a prismatic float with 
such beam loeiding that the chines do not Immerse during Impact. 

For waves that give the severe design condition of full-length 
impact, conventional heam loadings are smiLl enough to cause the 
maxlmimi force to occur at drafts sufficiently small to make the 
effects of finite width and china flare secondary. Reference 1 
Indicates that for a conventional float neglect of the pulled-up 
how la Justified \Jhen the trim is 3° or greater. Althou^ for 
high- trim landings ' initial contact hy the afterbody may sub- 
stantially change the trim before the main forebody impact^ the 
neglect of afterbody loads is Justified because, dmtLng the'main 
impausi^ the shielding of the afterbody by the forebody due to 
depth of the step aM to keel angle is sxuih as to minimize the ‘ 
importance of afterbody loads. - 


Structural .. 

A simplified representation of primary elasticity of an 
airframe is shown in figure 1. A rigid lower mass mj, is considered 
to be connected by a massless spring to a rigid upper mass m 0 . 

In deteimlning the fundamental bending of airplane wings part of 
the wing mass must be Included in and part of the wing lift 

should be applied to m^. In the present paper the gravity force 
on each mass is assumed to be balanced by wing lift. 

The problem of determining the properties of the two-mass 
system so that it is representative of the primary elastic action 
of the airplane is rather .simple if it is assumed that during 
the impact the structure deflects with the shape of its fundamental 
mode of vibration. The requirements ai’e: 
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(1) The total mass of the simplified system must equal the 
total mass of the airplane In order that the proper nodal or 
center-of -gravity accelerations can he obtained. 

(2) The energy of vibration for the same amplitude of the 

hull and lower mass (relative to the nodal pointy must \>e the same 
for the two -mas 8 system as for the considered mode of the air- 
plane structure. ’ 

• (3) The mtural frequency of the two -mass system must be the 
same as the- frequency of the considered mode of airplane vibration. 

Equations \dilch permit determination ’of the masses atffl. Sprlrl^ 
constant of the simplified system so that It meets these require- 
ments are given in appendix B . These equations and the’ foregoing 
requirements are appllcablei for both landplanes and seaplanes. 

in the present paper the represented structural mode is considered 
to be devoid of vibration prior to the i Instant of impact. Thus, 
the ccmputatiOTis may represent either a first impact or a subsequent 
impact resulting from a bounce sufficiently high to cause aero- 
dynamic and structural damping to stop the vibration during the 
time the seaplane is In the air. The present paper does not give 
a representation of successive impacts,; such as might occur In 
seaway, which lead to accumulative or resonant effects. Available 
fll^t data indicate that a single heavy impact, such as that 
considered herein, is the primary cause; of struct\iral failures. 

The response of the two-mass system Is obtained In connection 
with the calculation of the time history of the hydrodynamic 
force, and from this result the complete. response of the repre- 
sented mode can be obtained by the slmpie procedure given in 
appendix B and demonstrated in the .section entitled "COMPARISON 
WITH EXPERIMENT',." The response of other modes to the force 
computed on the basis of the fundamental mode can be separately 
determined and superposed (reference 3 ) . In order to minimize, 
the complexity of the solu'tloh, however, the present investigation 
does not provide for taking into account the effect of the other 
modes on the hsrdrodyhamlc force . Although the other modes may. 
have a substantial effect on the local loads in the structure, 
the effect of these modes on, the hydrodynamic force is considered 
to be secondary as compared with the effect of the ftmdamental 
mode. 


If a large number of solutions for the two-mass system have 
been made in order to determine the effect of the fundamental 
'mode of different wings on the hydrodynamic force, the response 
of modes other than the ftindainental can be approximated from 
the use of a solution for a mass ratio and ratio of the time 
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period of the in^iact force to the natural period of the tvomass 
system representative of the considered mode. An example of such 
use to approximate the effects of wing torsion is given in the 
present investigation in a comparison of results of computations 
\fith experimental results • 

If the response of more then one mode is considered, the • 
structural and aerodynamic damping, which are not considered in 
the present paper, are important factors in determining the extent 
to which the maximum response of the different modes should he 
superposed vrithout regard to phase relationship. It is expected 
that -the effect of the damping he most important for the 

higher modes and that a result leading to conservative design will 
he obtained if damping is not ccaasidered and the maximums of the 
first two or three mode's are superposed without regard to phase 
relationship. 


RESULTS 


Solutions of Ihe equations in appendix A were .made for wide 
ranges of the pertinent variables . Time histories of the calcvilated 
nodal acceleration, or .hydrodynamic force in terms of the weight, 
are given in figures 2, 3/ aiid ^ for ratios of the sprung mass 
to the lower mass equal to O.25, O.60, 1.00, and 1.3^, respectively. 
Each figure is three -dimens ionalj the third dimension is '^1/% 

which is a ratio of the period of natural vibration to the speed of 
the intact. In representing the period of natviral vibration, t^ 
is taken as the time required for one -fourth of a cycle. The speed 
of the impact is represented by B^ing t^ equal to the time between 

initial contact' and maximum- acceleration for a rigid structure. If 
the time to reach maximum force for the elastic structure should be , 
used in defining tj,, discontln\iities in the time to reach maximum 

force would cause discontinuities in 1lie time -ratio scales of the 
plots. (See figs. 2 to 5») . ■ 

An expression for may b© obtained from the relation 
tn = ^fji and equation (A4) in appendix A. The equation for t^^ 
is as follows: 


H = 


Ct 


’©max 


(-)' 

\PS/ 


1/3 


V 


( 1 ) 
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vhere a 

time coefficient at instant of maximim ^celeration for 
%ax rigid tedy (O.678 for P = 22.5°, 7 = and t s 3°) 

The expressions for tj^ and tj^ may he used to determine 

that 


^ ■ ac— (^) 

%ax 

The oscillatory acceleration Is given in figures 6 to 9. 

These figures are the same aS the figures giving the nodal 
acceleration (figs. 2 to 5), except that the acceleration plotted 
is the difference between the hull acceleration and the nodal 
acceleration. The time hiptory ef the hull acceleration can. he 
obtained by summing the- twe pictg a partlcu].ar mass ratio . 

On the basis of the assumption t^t. during impact the airframe 
structure deflects in a particular mode, a time history of the 
acceleration of any point in the structiire is obtained f ro^ the 
results for the equivalent two-mass system by the followi.ng 
procedure: 

(1) From the deflection curve of the represented mode obtain 
the ratio of the deflection of the point of interest to the. 
deflection- of the hull. Both deflections are taken relatiye to 
the nodal point. 

(2) Multiply this ratio by the Q§clllatory acceleration 
given either by figure 7, 8, or 9 or by interpolation between 
.these figures- for the. mass ratio of the equivalent two-mass system. 

(3) Add result to -tee nodal ac<ielerati on given by figures 2 

to 5* 


Time histories of the acceleration given in figures 2 to 9 
are on a npndimensional basis. The, nondimens ioiial Coefficient*, 
tdilch contain velocity, weight, fluid density, and acceleration 
of gravity, were used in reference 1 in a comparison of theoretical 
data with impact data for a float having an angle of dead rise 

of 22| . 

Application to other angles of dead rise .- The function of 
the angle of dead rise can also be included in the nondlmensional 
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coefficients, but In the present InTestigation this function Js 
Isolated and treated as a factor for correcting the results 

presented for angle of dead rise of 22— to other, dead-rise angles 1 

The pertinent relationships between results for different angles 
of dead rise may be expressed as follows: 


t « 




1 

f(P) f(A). 

(3) 

gO) f(A)j^/3 



where 

f(P) function representing variation of, virtual mass for two- 
dimensional flow with angle of dead rise 


f(A) aspect-ratio factor including effect of angle of dead rise on 
aspect ratio; ratio of virtual mass for three-dimensional 
flow to virtual Biass for two-dimensional flow 

np acceleration at any point, either oscillatory, nodal, or 
total 


The shape of the force and acceleration curves for a given 
value -of tp/tj is independent of angle of dead rise, but the 
effect of aixgle of dead rise on t^, as given by relation ( 3 ), 
does enter into the determination of the value of tp/tj for a 
particular solution. After the value of tp/tj for a particul^ 
Solution has been determined, the acceleration and time values 

• • lO' 

for an angle of dead rise of 22| are proportioned by means of 

relations (3) and (4) to the corresponding values for the angle 
of dead rise used in determining tp/tj in- ordor to obtain the 
proper acceleration history. ' 

Altho*ugh adequate impact data have not been available for 

checking the theoretical equations for angle's of dead rise other 
\0 

than 22^ , the theory is equally applicable to planing floats. 

Study of planing data has shown that the functions of angles of 
dead rise used in equations herein are approximately correct for 
angles of dead rise ranging from 15 ° to 30 °* The functions are: 


( 5 ) 
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f (A) = 1 - -iaJLJL (6) 

2 tan P 

Until Improved fmctione of Eoigles of dead rise are Obtained, 
f\nictions (5) aiid (6) shoiild be substituted in equations (1) 
and (2) and in relation (4) to correct for angles of dead rise 
within the range from 15° to 30° • Bough approximation can be 
obtained by use of functions (5) and (6) for angles of dead rise 
greater than 30° but not for angles of dead rise much less than 15°. 
In reference 2 there is a discussion of the inadequacy of f(A) 
for aspect ratios which normally occur for small dead-rise angles. 
For very large angles of dead rise and moderate velocity the 
static forces, which are not considered in the present inveatiga-_ 
tion, become of greater importance . 

Approximate correction to other flight paths and trims . - 
The flight-path and trim condition^ of primary interest to, the 
designer tend to be independent of the variables in seaplane 
design and constitute a limited range within which approximate- ' > 
correction of results to different angles of flight path and 
trims can be made without necessity for repeated time-history 
solutions. For any particular combination of horizontal speed, 
rate of descent, and trim, the most severe impact load for most 
of the structure occurs when the seaway is such that the keel 
contacts a wave slope approximately parallel to it. The effective 
angle of flight path and trim for such an impact are defined 
relative to the Inclined wave slope; therefore, the trim which 
gives maximum force is zero. The largest flight-path angle, 
relative to the keel and to the critical wave slope, is also /. ■: 
associated -vrith the most severe force. The value of the largest ' 
fli^t-path angle is not so definite as the critical trim but 
tends to be independent of , variations in size and wing loading. 

The velocity of the wave should be considered in determining 
the contact speed and fli^t-path angle. 


The equations and method- of solution given in appendix A 
permit solution for different flight*path angles and trims; 
however, approximate correction of the results in figures 2 to 9 
to other positive contact angles can be made by assuming that 
the proportionate effect of the structural elasticity on the 
hydrodynamic force is solely dependent on the ratio tn/tj^. 


Curves, given in reference 1 show values of Cj 


^max 


for different 


flight-path angles and trims, which may be substibtxted In 
equation (2) to obtain the value of, tn/tj for different contact 
angles. In making the approximate correction, the solution 
presented herein for the obtained value of the ratio t^/t. should 
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"be \u 3 ed to approximate the shape of the curve giving the desired 
time'histoj^. The load or acceleration scale shotild he corrected 
to the different contact angles hy multiplying the present result 
hy the ratio, determined from curves given in both references 1 . 
and 2, of the load-factor coefficient for the different contact 
condition to the load-factor coefficient for the conditions of 3° 
trim and l4*^ flight -path angle considered herein. Correction of 
the time scale Involves a similar procedtira in vhich the time 
coefficient is used rather than the load -factor coefficient. 

Since the force curves for a rigid body are approximately the 
same shape for different angles of flight path and trims (reference l), 
the approximate method of correcting to different angles of flight 
path and trims would be -almost correct if- the structural elasticity 
did not affect the hydrodynamic -force ' curve « The percentage change 
in the force on the float due to elasticity is a fumction of the 
percentage change in draft caused by the elastic compression for a 
given center-of -gravity position. An indication of the validity' of 
the approximate coi-rection- is obtained by studying the extent to which 
the ratio of the spring deflection to the draft is constant for impacts 
of the same values of tj^/t^ at different angles of flight path end 

trims, 3?rom the expressions for t^/t^, C^, C^, and with the 

spring deflection assumed to be proportional to the hydi’odynamic force, 
this ratio may be represented by the expression. CiC^/C^, in which: 

all values are for the instant of maximum acceleration,. The variation 
of this expression with flight path for an impact of a rigid float 
• at 3° and 12° trim is given in figure 10,. Values, of and 0,^ 

used in obtaining this figure are given in reference Ij values of 
were obtained in conjunction with the data of reference 1 but have 
not been published. 

In the present paper the numerical values of CiC^/C^ have no 

significance and they are of interest oray because of the extent to 
which they are constant. Figure 10 indicates tliat for large f light - 
path angles and small trims the ratio is approximately Constant. 

The deviation from a constant value of this ratio is due to planing 
forces which exist in an obliq^ue ■ impact and become more. Important for. 
low fli^t-patli angles and high trims. The conditions of large flight- 
path angle and small trim previously adjudged to be of' primary Interest 
to -Uie designer constitute the ranges in which the deflection ratio is 
fairly constant and the approximate method of correction should give 
a fair degree of accuracy . The present solutions are considered to 
be for conditions suited to correction of the results to other conditions 
of greatest practical interest^ they represent a moderately severe 
combination of flight-path angle, wave alope, and trim chosen to 
facilitate correlation of the theory with an impact which resulted in 
substantial damage to a well-instrumented flying boat during flight 
tests • 
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The equations presented herein are not valid for zero trim; 
an assumption that the float is prismatic gives solutions of 
infinite wetted length and infinite force for zero trim. Solution 
for 3 ° trim and a prismatic float is much simpler than a correct 
solxxtion for 0 ° trim hecaxise necessity for consideration of 
bow shape is eliminated . The solution for 3° trim may be talcen 
as an approximation of the critical design load or, as illustrated 
in the following section in a comparison, of theory with experiment, 
as an empirical factor, vrhlch includes bow effects, may be used to 
convert values of acceleration and time for 3 ® trim to values for 
0 ° trim. 


COMPAEISON WITH EXIERIMENT 


Impact “basic data have not been obtained for suitable models, 
and most flight landing data have been Inadequate for tlie present 
study. Tlie only data which appeared suitable for this comparison 
are those which were obtained with a four -engine flying boat, the 
data for which have not. been published. Data v;ere obtained for a 
large number of' test landings, but only one of the impacts is 
very well suited to the present analysis This impact gave loads: 
sufficiently high to cause large effects of elasticity of the 
wings The impact occurred against the flank of a sizable wave 
(4 ft) end thus facilitated the use of results’ based on a pl^ar 
water surface. The present comparison of theory vrith experiment 
will be restricted to this impact since other impacts- involved 
more complex contact with seaway and gave less force. 

The horizontal speed, rate of descent, and trim were recorded. 
A large number of pressure instruments distributed in the hull 
permitted deteiauination. of the water surface relative to the hixll . 
Data recorded by these Instruments indicated that the wave slope in 
contact with the hull was approximately planar, that the trim 
relative to the wave slope was 0*^, that the result^f velocity, 
considering the speed of the mve, was 85 feet per second, and that 
the flight -path angle relative to the wve slope v/as l4°. 

Structural data available for ’the test f]ying boat are not 
adequate for the present analysis. Since this flying boat has the 
same number of engines and approximately the same gross weiglit and 
horsepower as a landplane for which a large amount of structural 
data is available, assumption. is made that the flying -boat wing has 
the same mode shape and mass' distribution as the wing of this land- 
plane. Use in equation (Bo) in appendix B of data for the landplane 
given, in reference 3 leads to the following mass ratio of the tvo- 
mass system representing ihe fundamental mode; 
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Based on study of the accelerations at several points in the test 
flying hoat during periods of relatively free vihration in which the 
fundamental wing tending mode appeared to te predominant, a natural 
frequency of -3.6 cycles per second was selected for use in the 
present example. The fundamental mode frequency of the wing of the 
landplane 'is 3.4 cycles pel"* second. (See reference 3.) ' 

Since the mass ratio is equal to 0.25, figures 2 and 6 are 
used to approximate the action of the fundamental mode.. Further, 
sutstltution of values for conditions for this impact in equation (2) 

results in use of the specific time history given for = 1.2. 

■ti 

Sutstltution of the contact conditions in the load and time coef- 
ficients fixes the load and time scaiJ-ee in an atsolute sense. 
Gorrcctlon from 3° trim to 0° trim may te made ty assuming that 
the shape of the time history is appz'oximately the same for toth . 
conditions. Empirical correction of the curves from 3® trim to 0° 
trim may te seen from reference 2 to require a 10-percent reduction 
in the acceleration. An analysis of data attained at the 
Langley impact tasln. for Impact at 0° trim Indicates that correc- 
tion of the time scale from 3° trim to 0° trim requires a 10-percent 
reduction in the time values. 

Evdrodynamlc force .- The nodal-point acceleration nj-, attained 
ty the procedure discussed herein, represents the hydrodynamic force, 
in multiples of the weight, applied to the flying toat. ' Since the 
experimental data do not provide meastirement of the hydrodynamic 
force as such, direct comparison of the theoretical forcer curve 
results with experimental restfLts is not permitted. Instead, a 
comparison of the theoretical response of the structure with the 
experimental response is necessary, and, if the agreement is 
adeqtiate, it may be concluded Ihat both the hydrodynamic and the 
structural actions are adequately represented. 

Before a study of the response of the structure is made, the 
theoretical effect of the response on the hydrodynamic force should 
be. observed. This' observation is made by comparing' the force curve 
obtained for a mass ratio of 0.25 with the force curve for a rigid 

structure. Both curves are included in figure 11. The curve for 

/ mq \ 

the case of a rigid structure { ~ a 0 ] was obtained from reference 1 

for 3° trim and corrected to the conditions of the present example 
as previously indicated. Comparison of the cxurves for mass ratios 
of 0 and 0.25 shows that in the presen;b example the theoretical 
effect of the structural elasticity on the maximum hydrodjaiamic 
force is to reduce it I5 percent. For the hypothetical condition 
of a concentrated wing mass located at a point in each semlspan 
of a massless wing structure, the conditions of the present 
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example vovdd give theoretical reduction in the maximum hydrodynamic 
force due to structural elasticity of 44 percent. This result is 
indicated by coraparison of the, maximumi of the curve in figure 11 for 
a mass ratio of I.36 with the maximxmi of the curve for mass ratio 
of 0. 


The curves in figure 11 show only reduction of the hydro- 
dynamic force becaxise of structural elasticity. It shoiild not 
be concluded, hovever; that the effect is al^^ays in this direction. 
Figures 2 to 6 show that in some cases the hydro^mamic force is 
Increased; the maximum Increase *vdiich was calculated was of the 
order of 12 percent. 

Hull acceleration .- By ccanbinlng the results for the tvro-mass 
case with ai^roxiraation of the pitching action on the basis of a 
rigid structvire, the following e^viatipn may be obtained for the 
hvdl acceleration at different longitudinal stations: 




= n^f + n^ 


(7) 


where 


hull acceleration 




6 distance from station to center of gravity 
k radius of gyration (12 ft in present example) 

I distance from resultant hydrodynamic force to center of gravity 

The forebody length of the flying boat is 31.75 feet. Because 
of bow effects, a length of 25 feet is assumed to have a rectangular 
loading for this zero-trim impact and the resultant force is 
located 12. 5 feet forward of the step which leads to a value of 
1 eq,ual to 8 feet. Use of the foregoing procedure to calculate 
time histories of the acceleration for t\ro stations in the hull 
■at \diich accelerometers \ 7 ere located gives the curves in figure 12. 
The ma ximum accelerations recorded at these stations are also shown; 
agreement with the computed maximum acceleration is good. The 
f\oll experimental time history is not included becavise the film 
speed was not great enou^ to permit accurate determination of the 
shape of -toe time history. : This factor, together with some 
uncertainty in dexining the exact instant of contact, prevents exact 
check of the time to reach maximum acceleration; thus, the experi- 
mental points in figure 12 are located at the theoretical time of 
maxlmimx acceleration. 
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Although the agreement of calculation with experiment in 
figure 12 ie good, this agreement can' he interpreted as confirmation 
of the elastic action of the structure. only to the extent that 
disagreement of hull accelerations computed on the basis of a rigid 
stnacture can he shown. Values of maximum accelerations computed 
for a rigid structure are included in figure 12. The disagreement 
with experiment is greater than for computations in which the 
elasticity of the structure is considered. The difference is small, 
however, heca,xise the effect of elasticity in increasing the hvill 
acceleration for a given force ie largely offset hy the effect of 
elasticity in reducing the hydrodynamic force. Since the proh’ahle 
accuracy of available data must he considered, experimental proof of 
the. theoretical effect of the elasticity on the h^^ll acceleration 
Cjanno't he claimed. Support for a conclusion in' this respect, .however, 
can he obtained by comparing the theoretical response of the Trtng 
with the experimental response. 3‘f the wing responds as assumed, 
the basic equations require that the hydrodynamic force and hvill 
acceleratioiis he as calculated. 

Elastic axis .- Use of the nodal and oscillatory accelerations 
of the representative two -mass system to predict accelerations along 
the elastic axis of the wing requires consideration of the fact that 

the elastic axis of the wing is not at the center of gravj.ty of the 

flying boat. An approximate correction may he obtained hy multiplying 
the results for the twormass case hy the factor f which is used in 

equation (7). This correction is not entirely consistent with that 

given hy equation (7), hut each approximation is considered more 
accvirate for its particular case. An improvement to the present 
correction T-rhich would change , the results 2 to 3 percent- might he made, 
hut the complication is not considered to he warranted. 

In the present example the value, of f is 0.86. Application 
of this factor, of equation .(B7) of appendix B, and of pertinent 
structural data given in reference 3 gives the cterves in figure 13 
in representation of the acceleration time histories of the elastic . 
axis of the wing for the hull. Inboard -engine, nodal, outboard -engine, 
and tip stations. 

} • ' ■ 

Wing torsion .- The torsion of. the wing, during, Imj^ct.may have 
substantial effect on the acceleration of engine and lia, cells, . 
masses forward of the wing. Use of the procedure which gave the 
acceleration time histories in figure I3. to calcvilate the accelera- 
tion at the engine gives a maximian acceleration, of 3g at the 
inboard engines and a maximum acceleration of • 5. 6g, at •the 
outboard engines . The relative magnitude of these values is in 
strong disagreement with tlae fact that dtiring this Impact partial 
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failure of the inboard-engine mounts occurred, but the outboard- 
engine mounts were not damaged. 

Since the structural data used in the preceding computation 
are for the actual fundamental mode, their use in' in^iact calcula- 
tions involves assvmption that the coupling between the torsion 
and bending corresponds to the coupling which exists in natural 
vibration. Actually, the torsional deflections in impact are 
determined not merely by the bending deflections, but also by the 
large nodal acceleration, which does not exist in natural 
vibration. 

In the present example the results which have been calculated 
for impact'. of a two-mass elastic system will be used to predict the 
response of the engines. The procedure for doing this is to select 
the proper solution and then to adjust the acceleration and time 
scales of the two-mass solution so that the maximum acceleration of 
the nodal point corresponds to the maximum acceleration of the elastic 
axis at the engine station. The acceleration of the upper mass of 
the two-mass system then represents the response .of the engines; 
however, because of the eccentricity of the impact, an increment 
must be added. 

Data obtained frcm Guggenheim Aeronautical Laboratory, California 
Institute of Technology include the torsional deflection of a station 
inboard of the outboard engine for a given moment applied at the wing 
tip. For a flying boat the absence of cut-oUt for the landing gear 
tends to give a stiff er wing; therefore, in the present example; 
deflection measured on the wing of the landplane slightly inboard of 
the outboard engine is considered applicable to the oufboard-engine 
station. Relative, deflections between the inboard And outboard 
engines are estimated as follows: 

(1) Torsional deflection at inner engine equal to 1 imit. due 
to each engine, or 2 units total 

(2) Torsional deflection between inboard and outboard engines 
equal to 2 units because of greater distance of flexure, boosted 
to 3 units because of increased flexibility of structure 

On the basis of data from California Institute of Technology =, 
the average of the static moments of the inboard and , outboard engines 
is taken as 22,000 foot— pounds. After determination of the torsional 
deflection at the outboard engine for this average moment applied at 
the wing tip, multiplication by the ratio I/2 gives a value for the 
static deflection of the inboard engine and multiplication by the 
ratio 5A a value for the static deflection of the outboard 

engine. ' • 
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An approximation to the calculation of the response of the 
engines to the total acceleration of the elastic axis is to . 
neglect dynamic Interaction between the engines and treat each 
engine as a single-mass oscillator having a natural . frequency 
determined by its static deflection. Such a procedure gives 
values of 7*6 cycles per second for the Inboard engine and 
5.1 cycles per second for the outboard engine. Use of those 
frequencies to dotenalne tj^ and division of tn by t^ (tj 
equal In this case to the time to reach maj?:imum acceleration of 
the curves in figure 13 for the station in question) leads to 
values of- tji/tj for the Inboai^i and outboard engines. The 
next step is to select a mass ratio which for these values of 
tn/tj! has a shape of the nodal acceleration time history, 
which approximates the shape of the acceleration time history 
for the elastic axis at the station In question. In the present 
example the mass ratio of 0.25 Is used._ Scale factors for both 
the load and time scales are determined so that the maximum, 
nodal acceleration for the two-mass solution will agree with the 
maximum acceleration and time to reach maximum acceleration of 
the elastic axis at the station in question. After these factors 
are applied to both the nodal and oscillatory curves for the 
selected mass ratio and time-period ratio, use of the results 
and equations (B7), (B5), and (B2) in appendix B to calculate the 
acceleration of the sprung mass of the two-mass system gives 
accelerations of the engines. Approximation and superposition 
of the pitching action on the basis of a rigid structure leads to 
the solid-line emves given in figure l4 for the accelerations of 
the engines. 

Also included in figure l4 is the design ultimate acceleration 
for the engine mounts. Comparison of the calculated engine 
accelerations with this value shows agreement of the calculation 
with the fact that partial failure of the inboard -engine mounts 
occurred but the outboard -engine moiants were not damaged. 

During the impact an accelerometer was located at the 
outboard- engine station intermediate between the elastic axis 
and the engine. A calculated time history of the acceleration 
at the accelerometer location is given in figure lUj this time 
history is based on, linear interpolation between the computed 
accelerations at the engine and the elastic axis at this station 
in .accordance with the proportionate distances involved. The 
figure also includes the recorded maxlmxmi acceleration at this 
point and shows good agreement of the computed acceleration 
therewith. 

Acceleration time histories for the engines, computed on the 
basis of a rigid structure, are included, in figure 14. The maximum 
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acceleratlone computed on this basis do not agree vlth the 
structural failures which occurred. Furthermore, the fact that 
the curves computed on the' basis of a rigid structure reach a 
maximum at the same instant of time is in strong disagreement 
with experiment, which in this respect is in approximate agree* 
inent with the computations for an elastic structure. 

Difference between the ratio t^/ti for the inboard -engine 
and outboard- engine stations is primarily responsible for 
difference in the calculated response of the engines. Most of 
the difference in this ratio for the two stations is not due to 
difference in the sprung-engine frequency but is due to the 
greatly different time to reach maximum acceleration tj of 
the elastic axis. Agreement with experiment of the response 
calculated at these., stations provides indirect confirmation of 
the acceleration time histories predicted for the elastic axis 
at these locations by the normal -mode method. 

For the impact, experimental data are not available for 
checking the tip acceleration, which is predicted on the basis 
of the normal -mode method, but the initial downvai'd acceleration 
and the 12 g, maximum acceleration shown in figure I3 for this' > 
station are in general agreement with results recorded in severe 
Impacts of other airplanes. Agreement of the computed hull 
■acceleration with experiment has already been shown; in an 
indirect manner all the cur'/’es in figiure 13 exhibit satisfactory, 
agreement with available experimental data. Although the response 
of the engines is different from the response assumed in 
calculating these curves. It appears that in practical use 'the 
two-mass solutions given herein can be interpreted on the basis 
of the normal -mode method to obtain both the response of the 
elastic axis and the hydrodynamic force. Further, it appears 
that in practical problems modification of these, two results is 
not required when accelerations of the engine different from the . 
accelerations predicted by the normal -mode method are determined; 
the merit of this statement should be independent of whether such 
modification is made by the method used herein or by another method. 

In the foregoing comparison agreement of calculation with 
experiment is obtained without consideration of the response of 
modes higher than the fundamental mode . If, when more data are 
available^ it is shown that the response of the higher modes can 
be determined by treatment parallel to that given the fxindamentsd. 
mode, the two-mass solution given herein can be used to predict 
their response by selecting a solution, for a mass ratio and 
natural frequency representative of the higher mode, which has a 
nodal acceleration curve of approximately the same shape as the 
hydrpdynamlc-force cuirve. 
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CONCLUSIONS ■ ■ • ,r.,i 


. Theoretic^ solution .of hydrc4yn,amlc.. impact, ef a hvpfl;..?^ 

. connected, ty .a spring to an upper, imss, , '.the resxdt s'. of , o^ci^.ations 
for ;,vide ranges of mass ratio and. natur^’ frequency^; ;and .jthe-i'use 
of these results in a oomparison of theoretical data .,vith teat 
- .data for, a fli^t-test landing impact indicated,, that. :... .. T 

, 1. In flying-hoat impact the effect of the structural Vesponse 
•on the; hydrodynamic force might he suhstantial, the shaj^ of the 
force time history mi^t he conside.rahly changed, and^j the. maximum 
.hydrodjnami.c force might he either reduced or dncreased . , ^ 

2. The greatest'. reduction in hydrodynamic force occ.virre.d. •: 
for the condition , of large mass ratio and low .value of spring 
constant. , 

. 3 .. The normal -mode .method .was a. practical means for detemining 
the equivalent two-mass system which represented thp' major elastic ' 
action of the airframe, for predicting the effect of this action 
of .the, hydrodynamic force, and for approximating accelerations 
along the elastic axis of the wing. , , . ' 

h. The acceleration of. engines contained in nacelles forward 
of the wing could not he computed on' the basis of coupling between 
torsion and bending as in the fundamental mode, hut .simple .treat- 
ment .of their response to .the combined translation and o.sciliatlon 
.. .of the elastic axis gave agreements with experimental accelerations 
and gave an explanation of partial failure of the inhoard-engine 
mounts during impact. .1 ■ 

Lea^gley. Memorial Aeronautical Laboratory ■ 

National Advisory Committee for Aeronautics' 

Langley Field,. Va . , March 17, 19lj-7 . ■ 
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APPENDIX A 

MATHEMATICAL EQUATIONS AND METHOD OF SOLUTION 
Equations of Motion 


The following equations of motion, which were derived from 
equation (30) in reference. 2, are for fixed- trim impacts of a 
rigid prismatic float connected hy a massless spring to a rigid 
upper mass: 

The acceleration of lower mass normal to water svu-face in 
feet per second per second can he expressed hy the equation 


■h = + ^1 pos t)^ + 


ms 


Ayt, 




+ (®s + (3% • yLo^) 


+ 2K^A cos T r y ^3 dt + 3AK2_^ cbs^t f [ 
^0 ^0 ""0 


Ay^3 + idl 


(Al) 


The acceleration of sprung mass normal to water surface in 
feet per second per second is expressed hy the equation 


ys = ^ |(Ayj ;,3 + m[,)‘ yL'+ 3 Ayj^^(yL + cos t (A 2 ) 


The acceleration of nodal point normal to water surface in 
feet per second per second is given hy the equation 
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The spring constant in pounds per foot of deflection can he 
expressed hy the equation 


t^ff^,(ii^mg)fn^ . . 

' mi, + ns 


(Al^) 


where 




sln'^^T 
0 cos T 


A , 0.82 - iVfl - -^SUL\ 

\2p J y 2 tan sin t cos-^Ty 


y_ velocity Of lower mass normal to water surface, ft/sec 
•yj^ draft normal to water surface, ft 
0 angle of dead rise, radians 

T angle of trim, deg 


yj^^ initial velocity norroal to water surface, ft/sec 
xLq initial velocity parallel to water surface, ft/sec 
mg sprung mass, slugs 
lower mass, slugs 

fjj natural bending frequency, cycles/sec 
p mass density of water, slugs/ cu ft 

Computing Directions , 

A sample data sheet is given as table I. In this table the 
numbers in circles refer to rows; the circled numbers under the 
row headings refer to computed values to be used, for the computa- 
tions. In the first column the. time is equal to zero, in the 
second column the time is At, and in successive columns the time 
is 2 At, 3 At, and so forth. Each row is cmputed in sequence for 



NACA TO No. 1398 


21 


any glron column 'before any row Is computed for the next column, 
except for the first column. In the first column the assumed 
value of Is the given value of . All other values are 

zero except . tho values for rows and (0) which in this 

case are also Lower-case p as a subscript on a row 

number refers to the congjuted value in the Indicated row of the 
-preceding column. 


Each row heading Indicates the general operation to be 
performed. When the data sheet is set up numerical values should 
be substituted in the row headings for the expressions that are 
constant for a specific case. The constants for the sample 
computation, given in table II, are as follows: 


K . 

Ki 


0, deg . . . 
t, deg . . . 
ho> fVsec . 
ft/sec . 


....... 

cycles/sec 
p, slugs/cu ft . 
At, sec . . . . . 


10?66U 

U.2U317 

133-919 

22.5 

3 

20 .6673 
82.1585 • 
715.217 
525.778 
3 

1.938 

0.005 


For most solutions a time increment At of 0.005 second is 
satisfactory. (The value of At may be varied by considering the 
time for a given mass to reach a in;^Tri nrnm acceleration.) 

The number of significant figures to be used should be chosen 
on the basis of the computing equipment available and the accuracy 
desired. Comparison of, a solution computed with four significant 
figures with results which had been obtained with six significant 
figures gave a difference of about 1 percent in the maximum 
acceleration. 


As the computations proceed, the lower-^iass acceleration 
(row(|^), the sprung-mass^cceleration (row^), and the nodal- 
point acceleration (row^, all In g units, should be plotted 
against time (row©) in seconds. Each time a column is completed 
the new points should be added to the plot. This plot is the 
only brief method of checking on the accuracy of the computation 
until sufficient solutions are obtained to permit cross-plotting. 



22 


MCA..TN Hp.^1398 

If tho points/ do not' lie 034 a.siaooth ctrnre. an eri^ iW''been inaie 
In the conipulbations. If ';an jerror Is in^eiev©^ yaluo con^inted ' 
theroafltpr contaliis the error and therefore .j^eat care imet: 1)6, 
tajcen. , /'v - /,• " 

The nimber of columns required for the cosiqnitatl^s^ 
cn thp inm^oh of coiv^ required for rows ^ and ^ ; tbr;p.asa.:,- 
thelr ros^ctiYe maxlmuins. At least threo colutnns should lea.'rL^^ 
computed ^^yond the column In >diich row reaches its maxi mum 
ralue.' " ' 
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APPEOTIX B 


BELATIONS BETWEEN T^TO-MASS S^^TEM AND REPRESENTED 

STRUCTURAL MODE 


The sum of the masses and shovm in figure 1 must 

equal the gross mass of the represented airplane in order to 
obtain the proper nodal acceleration. For the hypothetical limit 
condition in which the wing mass is concentrated at a single 
point in each semispan, is the actual hull mass ajid iTig is 

the actual wing mass. In order to take into account the more 
complex natvire of the structural action for a particular mode 
the determination of the ratio of these masses is necessary so 
that the vibrational energjr of the siu^plified and represented 
systems are equal for the same vibrational aiiq)litudesof m^ and 
the actual hull or point of force application. 

On the basis of the' theory of vibrations, for which equations 
are included in reference 3> the vibration energy E of the two- 
mass system is given by the equation 


where 


«PL 

vibrational amplitude of 
system 


relative 

to nodal point of 

<Ps 

vibrational amplitude of 
system 

ms 

rela-tive 

to nodal point of 


o> natural frequency 

Since for the two-mass system the node is at the center of 
gravity 


Combining equation (Bl), equation (B2), and the fact that the 
total mass m is equal to the sum of the, masses ms and 


the 
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following equation can "be obtained for the vibrational energy E 
of the two-mass system: 


0,2 5 mg ^ 


(B 3 ) 


The vibrational energy of the represented mode is a function 
of the spanwlse mass distribution and mode shape. On the basis 
that the semispan of the airplane is divided into J sections or 
stations, the vibrational energy of the mode can be written as 
follows : 


= I + . . . + (B 4 ) 

where 

m. mass of Jth spanwlse section; value for semi span doubled to 
represent the entire span 

(p. deflection of mass at jth span^se section relative , to 
nodal point (h^ + m^) 

h* deflection of elastic axis at jth spanwlse deflection relative 
^ to nodal point 

X chordwise distance from elastic axis to effective mass center 

aj torsional deflection at jth spanwlse station 

Equality of the vibrational amplitude of the lower mass of 
the simplified system to the, vibrational amplitude of the- hull or 
fuselage of the flying boat or airplane relative to the nodal 
point of the represented mode is expressed by . the equation 

<Pl = ^ (®5) 

where 

<Pj^ deflection of hull or fuselage of flying boat or airplane 
relative to nodal point of represented mode 

The requirement of equal energy of the simplified and represented 
systems for the condition expressed by^ equation (B 5 ) gives combination 
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of equations (B 3 ), and (B 5 ) to ottaln the following equation 

for the mass ratio of the two-mass system: 


5s ^ 

mr p p p 




(b6) 


Computation of the natural frequency of ^ng modes has received 
a great deal of attention in connection with study of wing flutter 
and need not he treated herein* Incidental to calculation of the 
natural frequency, a mode shape is attained which, together with 
knowledge of the mass distrihution, permits use of equation (b6) . 

In cases in which the wing has .been constructed, the mode shape 
and natural frequency may be determined experimentally. Equatioo (Ak) 
in appendix A of the present paper permits computation of the spring 
constant which for a given mass ratio of the simplified system gives 
the required natural frequency. • 


After the accelerations of the two masses of the simplified 
system have been computed, equation (A 3 ) of appendix A fixes the 
magnitude of the nodal acceleration- The difference between the 
nodal acceleration and the hull, fuselage, or float acceleration 
can be taken as a measure of the oscillatory acceleration. On 
the basis that the structure deflects in the mode used in determining 
the equivalent two-mass system, the acceleration at any point is 
given by the equation 


Hp = nj + no ^ (B7) 

The foregoing equations, with consideration of rotatory 
Inertia and energy, may be applied to the case of a tip float 
attached to a flexible wing if the stiffness and mass distribution 
are Iknown and a manner of structural deflection is assmed. 


O 
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TABLE II.- SAMPLE CCMPOTATIOB 
[Conetante are from appendix aJ 


Column 


Rov 

Row heading 

1 

2 

3 

h 

5 

6 

7 

— j 

8 ; 

W 

Time At, sec ( 0 . 005 ) 

0 

0.005 

0.010 

0.015 

0.020 

0.025 

0.030 

0.035 

0 

©P 

20.6673 

20.6673 

20 . 61*20 

20 . 52 &T 

20.2688 

19.8290 

19.1964 

18.3859 

0 

0 . 005 ) ♦ ©p 

0 

3.103337 

0.206579 

3.309380 

0.411128 

0.511025 

3 .608174 

0.701700 

© 


0 

0.010679 

0.01*2675 

3.095716 

0.169026 

0.2611147 

0.369^6 

0.1*92383 

0 


0 

0 .001101* 

0.008816 

0.029613 

0 .069491 

0.1331*53 

0 .224949 

0.31*550^ 

® 


0 

0 .000111* 

0.001821 

0.009162 

0.028570 

0.068198 

0.136800 

0.242441 

0 

(0 ♦ U. 23736 )^ © X ‘* 01.757 

0 

2661.07 

10612.5 

’23586.5 

1*0782.0 

60767 .1* 

81602 .9 

101246 

© 

© X 33 .‘^ 79 S 

0 

0.003817 

0.060967 

0.30671*2 

0.956518 

2.28326 

4 .58030 

8.11608 


121 * 0.99 ^ - (20.6673 X ©)] 

0 

0 

-0.116653 

-0.781821* 

-2.75252 

-7 .02152 

-111. 8095 

-26.8740 


[(© ♦©?) °-°°«] ♦ ®P 

0 

0.000003 

0.000028 

0 .000124 

0.000373 

0.000879 

0 .001775 

0.003201 


113 ‘*. 92 x @ 

0 

0 .003U05 

0 

b 

. 

CD 

0.140730 

0.422190 

0.997595 

2.01448 

3.63288 

© 

[{© ♦ ®p) 0-«>25] . @p 

. 0 

0.000027 

0.000160 

0.000506 

0.001168 

0.002243 

0.003821 

0.005977 


[(© * ©p) o-®®®’] * ©p 

0 

0 

0 

0 .000002 

0,000006 

0 .000015 

0.000030 

0.000054 

@ 

7213.62 X (g) 

0 

0 1 

0 

0.014427 

0.043282 

0.108204 

0.216409 

0.389^*25 

® 

® ♦ ® ♦ ® ♦ © 

0 

0.007222 

-0 .023903 

-0.319925 

-0.133053 

-3.63246 

-7.S883I 

- 14.7356 

@ 

(150.533 X ©) *® 

0 

2662.16 

10608.9 

23538.3 

40581.7 

60220 .6 

8 o 415.5 

9902B.O 

@ 

(133.919 X @) + 525.776 

0 

525.921* 

526.957 

529.742 

535.082 

543.648 

555 .901 

572.91*6 

@ 

-@/@ 

0 

- 5 . 06 li 3 r 

-2O.I32U 

- 44.4335 

- 75.8420 

-110 .771 

-144 .658 

-I73.LI2 


(@ ♦ @p) 0.0025 

0 

-0.012655 

-0.062986 

-0 . 161415 

-0.300690 

- 0.466533 

-0 .638573 

-0 .794425 


© ♦ ®p 

0 

-0.012655 

- 0 . 075 ^ 1*1 

-0 .237056 

-0.53771*6 

-1 .00428 

- 1.64205 

-2.43728 

@ 

20.6673 + 0 

20 .6673 

20 .651*6 

20.5917 

20.4302 

20.1296 

19.6630 

19 .0245 

18.2300 



20.6673 

20 . 6 U 83 

20.5602 

20.31*95 

19.9793 

19.4297 

18.7052 

17.8328 


© * ? 

20.6673 

1 20 .6420 

20.52&7 

20.2688 

19.8290 

19.19<3l* 

18.3859 

17.4355 

@ 

[- 0.001390 (0 - 

0 

O.OOOOU7 

-0.000156 

i -0.00209: 

( -0.008696 

; -0 .02371*0 

> -0.05155s 

> -0.096297 

@ 

0/32.2 

0 

- 0 .157203 

-0.62523c 

) -1.37992 

-2.35531* 

- 3.44009 

-4 ,49248 

-5.37615 


1 (525.776 X ©) + (715.217 X ®) 
121 * 0.99 

0 

-0 . 06657 ! 

5 -O.26U98: 

J -0.58581*1 

» -1.00291 

-1.47116 

-1.93306 

- 2.33323 
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Fig. 1 
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Fig. 3 







Fig. 4 


NACA TN No. 1398 




NACA TN No. 1398 


Fig. 5 






Fig. 6 
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Fig. 7 







Fig. 10 


NACA TN No. 1398 



F/gure lO.-Vanai-jon of force- pens frof /on rot to at 
rr>axirr)um force w/fh flighf^pafh a/ig/e and tr/zoo. 




Fig. 11 



SUO4. ' 



Figure f I - Variation cf h'^dro dynamic fbrce with time* 
T-0° i 14'’ j \/=d& feet per second j W= 40,000 pounds j 

fi^=-3.6 cycles per second. 



Fig. 12 
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Accelerafion, At /sec 
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Fig. 13 



Time , sec 


Figure /3.~ i^ariofton of acceleraiions ai- stations along 
the r/mg elastic axis iv/th fume. 




-/o 


Fig. 14 
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Figure 1 4-. Fanaffon oi^ accelerations of engines irith time. 



